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ABSTRACT [In this work, we perform DFT + U periodic calculations to study geometrical and electronic 
structures and oxygen vacancy formation energies of Sm,CayCe1-x-yO2-5 systems (x = 0.0312, 0.0625, 0.125 and 
0.250; y = 0.0312, 0.0625, 0.125 and 0.250; 6 = 0.0312, 0.0625, 0.125, 0.250 and 0.50) with different oxygen 
vacancy and doping concentrations. The calculated results show that the V|-Sm**-V> structures where there is 
a position relationship of the face diagonal between Vi and V2 both nearest to Sm++ have the lowest energy 
configurations. The study on electronic structures of the Sm«CayCe1-x-yO2-5 systems finds that excess electrons 
arise from oxygen vacancies and are localized on f-level traps of their neighbor Ce, and Ca’ and Sm** 
co-doping effectively restrains the reduction of Ce“. In order to avoid the existence of Ce**, x and y must be 
both larger than 0.0625 as ô = 0.125 or 6 must be smaller than 0.125 as x = y = 0.0625. The Ce**/Ce** change 
ratio k has an obvious monotonous increase with increasing the vacancy oxygen concentration. The 
introduction of Sm°* decreases k. In addition, the doped Sm?“ can restrain the reduction of Ce** when the 
V1-Sm**-V> structure with a face diagonal position relationship in lower reduced atmosphere exists. It need be 
pointed out that the Smo.2sCeo.7501.5 system should be thought of as a Sm-doped Ce203 one. 
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1 INTRODUCTION 


CeO-based materials are important in solid oxide fuel cells (SOFC)?! due to their high ionic 


conductivity. Pure CeO2 is a mixed ionic and electronic conductor. When a neutral oxygen vacancy is built in 


CeO2, two excess electrons left behind will be localized on the flevel traps of two neighbor Ce, namely, 
formally Ce** is reduced to Ce**[ Cep, ]'**1. Therefore, to improve its ionic conductivity and decrease its 


electronic conductivity'”!>!, CeO2 is commonly doped with alkaline earth or rare earth metals such as Ca**, 
Gd”, Sm* and La% or their co-doping. 

Now, ionic and electronic conductivities of CeO% have been widely studied!”-'®!. Mogensen’s study on Ca 
or Gd single-doped CeO systems found that the conduction carrier was related with atmosphere!'*!. The 
conductivity measured in air is regarded as being ionic and the electronic conductivity contribution is 


negligible, and ionic conductivity (oi) and electronic conductivity (oe) are connected with oxygen partial 


pressure ( Po, ) in reduced atmosphere!*!. The highest oi of CeO2 doped with 20mol% Ca is 0.209 S/cm at H2 


atmosphere at 850 °C, larger than that in air!!?!, Cei..M,O25 (M = Gd** or Sm**) at x = 0.15, 0.2 and 
(Ceo.sSmo.2)1-0.0sCa0.0s02-5 compounds have the largest ionic conductivity], and the conductivity of the latter 
is 0.126 S/cm at 800 °C in air. In addition, the ionic radii and valence of dopants!!®°), oxygen vacancy and 
dopant concentrations”!!, and dopant-vacancy and vacancy-vacancy defect associations'*”! which decrease the 
number of mobile vacancies available, also affect the ionic conductivity of CeOz. 

Theoretically, the distributions of oxygen vacancies and dopants"), their formation energies'**!, geometric 
and electronic structures'?>! and oxygen ionic migration energies of CeO. systems!!! 25: 761 have been 
investigated at the atomic level. CeO2 systems with Sm** among rare earth metals and Ca?“ among alkaline 
earth metals have the smallest oxygen ionic migration energies!*>-*7! and their doping can effectively restrain 
the electronic conductivity. Theoretical”! and experimental!?* °°! studies showed that Ca?“ and Sm** co-doping 
could better improve the ionic conductivity. Effects of different oxygen vacancy and doping concentrations on 
the distributions of oxygen vacancies, dopants and Ce** of the CeO: systems are investigated. Murgida’s 
study®°! showed that the oxygen vacancy concentration affected the distribution of Ce?“ and excess electrons 
preferred to be localized in the cation sites such that the mean Ce** coordination number was maximized, and 
two vacancies were inclined to be second-neighboring. For doped CeO2.s, our and other studies found that 
dopants preferred to occupy the nearest neighbor (NN) and next-nearest neighbor (NNN) positions relative to 
the oxygen vacancy! 251, Independent with oxygen vacancy and doping concentrations, Ce** is distributed 
around the nearest neighbor oxygen vacancy°*71, 


However, a systematically theoretical study on Ce**/Ce** electronic properties of CeO2 with different 


oxygen vacancy and doping concentrations is absent. As is well known, Ca and Sm doping can effectively 
restrain the reduction of Ce** to Ce** and the doped CeO; has the largest ionic conductivity. Hence, in order to 
better understand the effects of different oxygen vacancy and Ca-, Sm-doping concentrations on the electronic 
structures (Ce**/Ce**) of CeO systems, in the work, we select various numbers of oxygen vacancies or Ca and 
Sm dopants in the same supercell to obtain a series of vacancy and doping concentrations. Ultimately, this 
work would plot transitional area of Ce** to Ce? and the Ce**/Ce** change ratio k under different oxygen 


vacancy and doping concentration conditions. 


2 MODELS AND COMPUTATIONAL DETAILS 


2.1 Models 

CeO; has a fluorite-type structure (Fm3m space group) with one formula unit per primitive unite cell!?>!. 
In this work, we considered a number of oxygen vacancies or dopants of Ca*+ and Sm** to produce various 
oxygen vacancy and doping concentrations in a 2x2x2 supercell, namely, Sm,CayCe1.x-,O2-5 systems with 
various x, y and 6 (x = 0.0312, 0.0625, 0.125 and 0.250; y = 0.0312, 0.0625, 0.125 and 0.250; 6 = 0.0312, 
0.0625, 0.125, 0.250 and 0.50) were considered, see Table 1. When the doping concentration is larger than 0.30, 
Sm,CayCe}.x,O2.5 systems were not doping but alloying ones. However, for convenience with description in 
the work, we would unify them to name as doping ones. 

Doped atoms are uniformly distributed in a 2x2x2 supercell to obtain systems with various doping 
concentrations. Hooper’s study?! on Sm-doped CeO: systems found that the dopant-vacancy interaction was a 
hybrid NN/NNN mixture distribution as the Sm** concentrations increased and more NN distributions were the 
most favored. Our previous studies found that the first oxygen vacancy (V1) was the nearest to the dopant 
(NN)? 25]. Here, the same structure model is used. For the CeO systems where the number of the oxygen 
vacancy nearest to one Sm** is more than one, these vacancies are chosen according to the rule of the NN 
distribution, namely, the second oxygen vacancy nearest to one Sm** is introduced and named as V2. Thus, Vi 
and Vz would have three distinct position relationships in the cube of eight coordinated O% for one Sm**, 
namely, the side one in Fig. la, the face diagonal one in Fig. 1b and the body diagonal one in Fig. 1c. 
Calculation results show that the structures with the face diagonal relationship between Vi and V2 are the most 


stable, consistent with the case of pure CeO» systems with double oxygen vacancies'®!, For Ca-doped CeO2 


systems, the models where Ca? has a NN or NNN distribution relationship with one oxygen vacancy are 


considered. 


(a) (b) (c) 
Fig. 1. Geometric structures of Smo.25Cao.25Ceo.sO02-5 or Smo.sCeo.sO2-5 systems with Vi and V2. Red spheres note O, 


white ones note Ce, and the pink ones note one Sm**. There is a similar notation in the following figures 


2.2 Computational details 

All calculations were performed by a Vienna ab initio simulation package (VASP)4!. Ce5s5p6s4f5d, 
O2s2p, Ca3s3p4s and Sm5s5p6s4f5d were treated as valence electrons. Structures were relaxed until forces on 
each ion were below 0.02 eV/A and the total energy was converged within 1x10“ eV. A plane-wave cut off 
energy of 400 eV, a 3x3x3 Monkhorst-Pack k point mesh and a Gaussian smearing parameter of 0.20 eV were 
used. 

The standard DFT formulation usually fails to describe strongly the correlated electrons due to a deficient 
treatment of electron correlation. This limitation can be corrected by using a DFT + U method, where the 
introduction of a Hubbard parameter U modifies the self-interaction error and enhances the description of the 
correlation effects% This methodology has been widely used in reduced CeO2 systems'®: * 32-351. Theoretical 
work showed that the U value for Ce should be larger than 5.0 eV>- 26, 341, In this work, we used the U-value of 
6.0 eV for Ce and U-value of 8.0 eV and J-value of 0.65 eV for Sm, consistent with the other work!?>-47!, The 
exchange-correlation effects were described with the Perdew Burke Ernzerhof (PBE) functional within the 
generalized gradient approximation (GGA)*!. The calculated crystal lattice constant from the GGA + U 
method is 5.48 A, in agreement with the experiment value of 5.41 Al8), 

The formation energies of oxygen vacancies Evo for CeO2.5 systems, Evo-sm for Sm,»CeO2.5 systems and 


ENo-Casm for Sm,CayCeO2.5 systems can be defined as 


) = E[X ]+2 E[O,]-E[Y] 


Ey, —(Sm,CaSm 


(1) 


where E[X] and E[Y] are total energies of pure or doped CeO2 systems with and without oxygen vacancies, 
respectively. E[O2] is the energy of one O2 molecule set in a 10A x 10A x 10A supercell, and p is the number 


of oxygen vacancies. 


3 RESULTS AND DISCUSSION 


3.1 Geometric structures 

Model structures with the lowest energies are obtained. The Vi-Sm*+-V> structures where there is a 
diagonal position relationship of the face between Vi and V>2 both nearest to Sm** have the lowest energy 
configurations, consistent with the case of pure CeO systems with double oxygen vacancies'*!. The 
introduction of Sm*+, Ca?* and oxygen vacancies into CeO2 can produce obvious geometric distortions, 
consistent with the other work’, Geometric structures of CeOj.96ss, Smo.0312Ceo.96880 1.9688, 
$m0.0312C 0.96880 1.9375, SmM0o,.0312Cao.0312C 0.93760 1.9688 and Smo.0312Cao.0312C 0.93760 1.9375 Systems as examples are 
displayed in Fig. 2. Geometric structures of Sm,CayCe1-x-yO2.5 with other oxygen vacancy and doping 
concentrations have a similar geometric distortion. 

It is well known that the oxygen vacancy is an area of effective positive potential, hence, the neighboring 
O> move toward the vacancy, and the neighboring Ce** move away from the vacancy. From Fig. 2a of the 
CeO} 19688 system, we can see that three of four O% near the oxygen vacancy move toward it and another O” 
moves away from it. From the following electronic structures (see detail discussion on electronic structures), it 
is known that this O” is bridged by two Ce?“ which have larger negative potential than that of Ce**. From Fig. 
2b of the Smo.0312Ceo.96880 1.9688 system, it is similar to the CeO1.9688 system in Fig. 2a and the difference is that 
one of two Ce** is Sm”, namely, the neighboring O% move toward the vacancy, and the neighboring Ce** and 
Smt move away from the vacancy. From Fig. 2c of the Smo.0312Ceo.96ssO1.9375 system, we can see that the 
movements of four O2 toward two vacancies and two O* toward Sm** resulted from the common attraction of 
Vi and V2 to O”, namely, an O? bridged by Sm** and Ce** is not repelled, different from the systems with one 
oxygen vacancy. From Fig. 2d of the Smo.0312Cao.0312Ceo.9376O 1.9688 system, four O2 near the oxygen vacancies 
can be driven toward the vacancy and the neighboring Ceft move away from the vacancy. The 


Sm0.0312Cao.0312Ce0.93760 1.9375 System in Fig. 2e is similar to the Smo.0312C€0.968801.9375 system in Fig. 2c. 
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(d) 


(e) 


Fig. 2. Optimized geometric structures of the CeO1.968s (a), Smo.0312C€0.968801.9688 (b), Smo.0312C 0.96880 1.9375 (€), 


Smo.0312Ca0.0312Ce0.937601.9688 (d) and Smo.0312Cao.0312Ceo.937601.9375 (€) systems. 


Arrow directions indicate moving ones of ions 


3.2 Oxygen vacancy formation energies 


Table 1. 


CeO.5 
ô Ev 
0.0312 2:57 
0.0625 2.56 
0.125 2.81 
0.250 3.66 


Oxygen Vacancy Formation Energies (Unit: eV) of Evo for CeO2-s, 


Evo-sm for SmxCe1-xO2.5 and Evo-casm for SmxCayCet-x-yO2-5 


Sm,Ce}.,O2.5 


0.0625 


0.125 


0.250 


Sm,Ca,Ce1-x-y02-5 


0.0625 


0.125 


0.250 


Evo-casm 


Oxygen vacancy formation energies of Evo for CeOz.5 systems, Evo-sm for Sm,Ce1..O25 systems and 
Evo-Casm for Sm,Ca,Ce1-x-yO2-5 systems are listed in Table 1. Evo monotonously increases with increasing 0, see 
Table 1 and Fig. 3. For Sm,Ce-.O25 and Sm:CayCe1-x-yO2-5 systems, we can find that Evo-sm with a certain x 
and Evo-casm With certain x and y are large as ô is large, and both with a certain ô are small as x or x and y are 
large, see Table 1 and Fig. 3. It need point out that the introduction of Ca and Sm makes the oxygen vacancy 


spontaneously form, similar to Fergus’s study on the Sm-doped CeO: systems'*!, 
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Fig. 3. Variation of Evo, Evo-sm, and Evo-Casm as 6 of corresponding CeO2.s, Sm,Ce1..O2-5 and Sm,CayCe1-x-yO2-5 systems 


3.3 Electronic structures 

When oxygen vacancies in CeOz are formed, the Ce4f state is split into two states: an empty Ce4fempty 
state and an occupied defect Ce**4fjun state at the range of O2p and Ce4fempty, consistent with our previous!?>! 
and the other work'?!, see Fig. 4. Total electronic densities of states (DOS), partial electronic densities of states 
(PDOS) and localization electronic densities of states (LDOS) from the defect state of Cet for various 
Sm,Ca,Ce1-x-yO2-5 systems with different x, y and ô are calculated, as shown in Figs. 4 and 5. 
3.3.1 Excess electron distribution 

From the PDOS of Ce4f state for CeO1.9688, CeO1.9375, CeO1.375 and CeO1.75 systems in Fig. 4a, we can 
see that a new peak appears at the range of —1.2~0 eV for CeO1.9688, CeO1.9375, CeO1.875 systems and —0.80~ 


0.60 eV for the CeO 175 system, respectively, which are fully occupied by Ce** electrons. Compared to the 
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PDOS of Ce4f state for the CeO1.75 system, the Fermi level of the other systems approximately shift up by 0.50 
eV, due to the decrease of oxygen vacancy concentration. Excess electrons arise from the oxygen vacancy and 
are localized on the f-level traps of its neighbor Ce, which can be visually recognized from the corresponding 
LDOS of the defect Ce** state in Fig. 4b-e. These are consistent with theoretical studies for CeO; 943, CeO1.9687, 


Ce0 9375 and CeO;.875'*3) and CeO2(111), (110) and (100) surfaces", 


DOS(states/eV) 


PDOS\(states/eV) 


(b) ô=0.250 


(d) ô=0.0625 (e) d=0.0312 


Fig. 4. DOS of the CeO2 systems and PDOS of Ce4f state for CeO1.9688, CeO1.9375, CeO1.875 and CeO1.75 systems (a) and 
the corresponding LDOS (b)-(e) of the defect peaks. The isosurface is shown in green and is set to 0.05 e/A. Here, one 


primitive cell of various CeO2.5 systems with a 2 x 2 x 2 supercell is shown 


From the PDOS of Ce4f state for Sm,Ce1-.O2.5 (0.0312<x<0.25, 0.0312<d<0.5) systems, see Fig. 5a, 


similar to CeO2.s systems mentioned above, a new peak appears in the range of O2p~Ce4fempty. From the 
corresponding LDOS of the defect Ce? state, see Fig. 5c-j, we can see that excess electrons arise from oxygen 
vacancies and are localized on f-level traps of their neighbor Ce, like the case of CeOz.5 systems!>: 26, 30, 321, 
For the Smo.0312Ceo.968g01.9375 system, see Fig. 5j, there are two oxygen vacancies in the 2x2x2 supercell, so 
they should induce four charge-compensation cations. However, the calculated result finds that there are three, 
maybe due to the existence of Vi-Sm**-V> structure with a face diagonal position relationship in lower reduced 
atmosphere and then doped Sm** can restrain the reduction of Ce**. 

In order to better restrain excess electrons, Sm,CayCe}-x-)O2.5 (0.0312<x<0.25, 0.0312<y<0.25, 
0.0312<d6<0.5 ) systems with various x, y and 6 are explored. From the DOS of Fig. 5b of Sm,CayCe1-x-yO2.5 
systems, we can see that, except for the Smo.0625Cao.0625Ceo.g750 1.875 system, Sm,CayCe1-x-yO2-5 systems have no 
Ce**. In other words, in order to avoid the existence of Ce*+, x and y must be respectively larger than 0.0625 as 
ô = 0.125 or 6 must be smaller than 0.125 as x = y = 0.0625. For a series of Sm,Ca,Ce1-x-yO2-5 systems, from 
their corresponding LDOS of the defect Ce** state in Fig. 5c-k, we can see that excess electrons arise from 


oxygen vacancies and are localized on f-level traps of their neighbor Ce, like the case of CeOo.9!> 26: 30. 321, 


sE s{(al)x=y=0.250, &0.250 ; 
af — 1=0.0312, 60.0312 | ‘ i S | 
Ki — 00312, 0.0625 : : 7 N | 
1 F(b1)x=)=0.250, O=0.50 
26 J 6 | 
1E 4 | 
; i CA P dhan A = x. 
oF i l Mel) x=y=0.125, 20.125 ry 
gh — 190.0625, 60.0625 : s [NA 
anyi — 30,0625, 620.125 S : | 
> a "idly x=y=0.125, 0.250 
9 7 Be VAS 
Fi | a) N 
$ 0 Ba i \ 
§ st a 1 Kel amy0 0625, 00625 VT | 
Zat = x=0.125, &0.125 Ni | 
3 si 120.125, 0.250 Q: 
A a a [| (£1) x=y=0.0625, 60.125 | 
Au IE N 
Ha i (g1) x=y=0.0312, 60.0312 | 
4h K0250, 60.250 : | | 
aE 120.250, &0.50 : | | 
ab Hal) x=y=0.0312, 320.0625 | 
1È ‘ | | 
2 | hn 
be 5 -20 15 -10 5 i - 5 
E(eV) 
(a) (b) 


D) a & 
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(i) x = 0.0625, ô = 0.125 (j) x = 0.0312, ô = 0.0625 (k) x = y = 0.0625, ô = 0.125 
Fig. 5. PDOS of Ce4f states for the Sm+Ce1-x0O2-5 (0.0312 <x <0.25, 0.0312 <8 <0.5) systems (a) 
and DOS of SmxCayCe1-x-yO2-5 (0.0312 <x <0.25, 0.0312 Sy <0.25, 0.0312 <8 <0.5) systems (b) 


and the corresponding LDOS (c)~(k) of the defect peaks. Oxygen vacancies are represented 


by Vi and V2. The isosurface is shown in green and set to 0.05 e/A 


3.3.2 Transitional area of Ce to Ce** and Ce**/Ce* change ratio k for Sm,CayCe1.x-yO2.5systems with 
different oxygen vacancy and doping concentrations 

In order to visually understand transitional area of Ce** to Ce** for the Sm,Ca,Ce1-x-yO2-5 systems with 
various oxygen vacancy and doping concentrations, their transitional areas of Ce** to Ce? for the Sm,Ce1-x02-5 


(Fig. 6a) and the Sm,Ca,Ce1-x-yO2-5 systems (Fig. 6b) are plotted. 
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Fig. 6. Transitional area of Ce** to Ce% and Ce**/Ce* change ratio k for Sm,Ce1.,O2.5 and 


Sm,CayCe1-x-yO2-5 systems with different oxygen vacancy and doping concentrations 


From Fig. 6a, we can see that no Ce** exists for the Sm,Ce1.,O2.5 systems with x = 0.167 and ô < 
0.0833 and Ce** exists for the Sm,Ce1-xO2-5 systems with x < 0.0312 and ô = 0.0312, x < 0.0625 and 6 
= 0.0625, x < 0.125 and ô = 0.125, x < 0.167 and ô = 0.166, and x < 0.250 and ô = 0.250. For 
the Sm,CayCe1-x-yO2-3 system as y # 0, the substitution of Ce** by Ca°™ makes two excess electrons and the 
substitution of Ce** by Smt makes one excess electron. Based on the case, for convenience of totally 
reflecting the effect of the doping concentration on the transition ratio of Ce** to Ce**, the doping effect of one 
Ca” is transformed to that of two Sm**, in which the ionic radius of doping Ca**+ and Sm** is omitted, and its 
corresponding transitional area of Ce** to Ce** is plotted in Fig. 6b. From Fig. 6b, we can see that no Ce** 
exists for the Sm,Ca,Ce1-x-,O2.5 systems with x = 0.0936 and ô < 0.0625, x = 0.334 and ô < 0.166, x 
= 0.375 and ô < 0.250, and x = 0.750 and ô < 0.50 and Ce** exists for the Sm,CayCe}.x-yO2-5 systems 
withx < 0.187 andd = 0.125. 

The Ce**/Ce** change ratio k in Sm,Ce}.,O2.5 systems with various oxygen vacancy and doping 
concentrations are studied, as shown in Fig. 6c. From Fig. 6c, we can see that k has obvious monotonous 
increase with increasing the vacancy concentration, and the introduction of Sm** reduces k, such as 
Smo.0312Ce0.968801.9375, Smo.0625Ce0.937501.875, Smo.12sCeo.87501.75, Smo.166Ceo.s3401.668 and Smo.25Ceo.7501.5 


systems, namely, the introduction of Sm** restrains the reduction of Ce** to Ce*. It need be pointed out that 
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Smo.25Ceo.7501.5 system can be thought of as Sm-doped Ce203, because of no unoccupied defect states between 


the occupied Ce** states and the unoccupied Ce4fempty, see Fig. 5a, in accordance with other theoretical 


workl51-531, 


4 CONCLUSION 


In this work, the influence of Sm** single-doping or Ca? and Sm** co-doping, oxygen vacancies and their 
concentrations on the geometric and electronic structures of Sm,Ca,Ce1-x-yO2-5 systems are studied. Results for 
the geometric structures show that the V;-Sm**-V> structures where there is a position relationship of diagonal 
of the face between Vi and V2 both nearest to Sm** are easily obtained, in agreement with the case of pure 
CeO systems. Study for electronic structures finds that the oxygen vacancies are contributed to excess 
electrons and these electrons are localized on f-level traps of its neighbor Ce for the Sm,Ca;Ce1-x-y02-5 systems. 
The Ce*+/Ce* change ratio k is related with the dopant and oxygen vacancy concentrations. In addition, the 
existence of V|-Sm**-V> structure with face diagonal position relationship in lower reduced atmosphere maybe 
makes the doped Sm** restrain the reduction of Ce**. It need be pointed out that Smo.25Ceo.7501.5 system can be 


thought of as Sm-doped Ce203. 
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The Ce*+/Ce* change ratio k has obvious monotonous increase with increasing the vacancy concentration. 
The introduction of Sm+* reduces k. It need be pointed out that the Smo.25Ceo.7sO1.5 system can be thought of as 


Sm-doped Ce203. 
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